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TECHNICAL MEMORANDUM NO, 1138 

ON THE APPLICATION CF THE ENERGY METHOD TO STABILITY PROBLEMS 1 

By Karl Marguerre ' - \ . 


• Since stability problems Have come into the field of vision of engi- 
neers, energy methods have proved to he one of the most powerful aids in 
mastering them. For finding the especially interesting critical loads 
special procedures have. evolved that depart somewhat from those cust omar y 
in the usual elasticity theory. A clarification of the connections seemed 
desirable, especially with regard to the postcritical region, for the 
treatment of which these special methods are hot suited as they are. 

The present investigation 2 discusses this question — complex (made 
important by shell construction in aircraft) especially in the classical 
example of the Euler strut, because in this case — since the basic fea- 
tures are not hidden by difficulties of a mathematical nature — the prob- 
lem is especially clear. 

The present treatment differs from that appearing in the Z.f.a.M.M. 
(1938) under the title "tiber die Behandlung von StabilitHtsproblemen mit 
Hilfe der energetischen Methode" in that, in order 1 to work out the basic 
ideas still more clearly, it dispenses with the investigation of behavior 
at "large" deflections and of the elastic foundation; in its place the 
present version gives an elaboration of the 6th section and (in its 7th 
and 8th secs.) a new example that shows the applicability of the general 
criterion to a stability problem that differs from that of Euler in 
many respects. 

5— n 1 1 1 " 1 ■ — — 

"Uber die Anwendung der energetischen Methods auf Stabilit'ats— 
problems, " Jahrb. 1938 DVL, pp. 252— 262. .. 

2 

In the paper investigations were continued at the instigation of 
Professor Trefftz (during his activity at the Deutschen Versuchsanstalt 
fur Luftfahrt). For a large part of the work (especially in secs. 4 and 
6) the author is very grateful to his colleague, R. Kappus, for his close 
collaboration. 

3 

See the , next to last paragraph of the Introduction, 
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SUMMARY 


In the two examples of the Euler strut and the slightly curved team 
under transverse load it was shown that the difference "between the sta- 
bility probl ems and the problems of linearised elasticity theory rests 
upon the fact that in the stability problems the expression for the energy 
of deformation contains terms of higher than the second order in the dis- 
placements. This idea makes it possible to establish the connection be- 
tween the energy method in the special form most used for stability in- 
vestigations and the principle of virtual displacements in its general 
elasticity - theoretical version; . besides, it permits the Investigation 
of elastic behavior beyond the critical deflection. 


INTRODUCTION 


K3rchoff , s uniqueness law states: An elastic body can assume one 

and only one- equilibrium configuration under a given (sufficiently small) 
external loading. In the formulation from the energy point of view: the 

potential II of the inner and outer forces has one and only one extremal 


511 = 0 


( 1 . 1 ) 


and the extremal is a minimum. 1 

The uniqueness law holds without restrictions in the realm of line- 
arized theory of elasticity, that is, as long as the stresses a, t can 
be expressed linearly in terms of the strains and the strains line- 

arly in terms of the displacements u, v, w. Then the function H is of, 
at most, the second degree in the displacements, and geometrical considera- 
tions show directly that a "parabola" of the second degree (positive def- 
inite quadratic form) can have one and only one minimum (or in mechanical 
terms, equilibrium position). The situation changes, however, when struc- 
tures are considered the behavior of which can no longer be expressed with 


1 Eor the derivation of the principle of virtual displacements (equa- 
tion (l.l)) for elastic equilibrium, see, for example,- reference 1, pp. 

70 ff. A careful foundation of the general theory of the behavior at the 
stability limit appears in reference 2, p. 160. (For further literature, 
see reference 1, pp. 277 ff., or reference 2). The investigations in 
section Stable and Unstable Equilibrium in particular majc® use of the 
Trefftz point of view. 
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sufficient accuracy "by the linearised strain-displacement equations. Such 
are, particulary, bodies for which One dimension is small compared to the 
others, structures in the form of shells, plates, or "bars. For example, 
a rod can, without exceeding the proportional limit, undergo bending de- 
flections several times greater than its thickness, and under these cir- 
cumstances the quadratic part of the (transverse) displacements in the 
strain displacement equations is no longer small compared to the linear 
part. Then the energy of deformation of 'the potential II becomes of 
higher than idle second degree in the displacements, and a parabola of • 
higher order can naturally have several extremals (equilibrium positions). 

The problem of the theory of stability is usually considered to be 
the determination- of that external load wider which several neighboring 
equilibrium configurations are possible. The reason for limiting in- 
vestigation to this "critical point" lies in the fact that the differential 
equations describing the elastic behavior in the postcritical region are, 
in general, no longer linear and an analytical treatment would therefore 
be difficult; while at the critical point itself the problem can still be 
"linearized". 1 

This purely practical viewpoint has, however, led to a certain (as 
is shown, from unfounded standpoint) systematic separation of the sta- 
bility problem from the other problems of the theory of elasticity, which 
finds its mathematical expression in a formulation of the principle of 
virtual displacements somewhat different from the usual one — has also 
for convei:ience led to the formulation of a special principle. (For ex- 
ample, see reference 3i). 

The principle of virtual displacements states that during a virtual 
(that- is, geometrically possible) displacement from an equilibrium posi- 
tion, ■ the energy of deformation taken up by the elastic body is equal to 
the work dene by the external forces. For the use of this principle in 
the theory. of elasticity it is convenient to express this fact in the 
following way: An equilibrium state is distinguished by the fact that 

for every virtual displacement from that state the potential of the inner 
and outer forces 


IT = A 3 + V 


Knowledge of the postcritical region was until now of secondary 
practical interest, because buckled structural, elements were considered 
unpermiS3ible. It has been only in recent years that in the shell con- 
struction of aircraft critical loads have "been permitted to be exceeded 
by large amounts unhesitatingly. 


4 


NACA TM No. 1138 


has. a stationary value: ■ 

5TI = 8^ + V) = 0 (1.2) 

Therein the potential of the Inner forces is given by the energy 

of deformation (inner /work), and the potential V of the .external forceB 
"by the negative product of the external forces considered constant, and 
the displacements of their points of application. In the region of appli- 
cability of the' proportionality law numerically V = ~2A a , where A a 
is the work done by the external forces as they increase from zero to 
their final values in passing through- only equilibrium states. The prin- 
ciple (1.2) can therefore be written conveniently, also in the form, (For ... 
example, see reference 3.) ' ' 

8(A i - 2A a ) « 0 ■ (1.3) 

As against this there is often used as a "minimal principle" in stability 
theory the condition 

StAi - A a ) * 0 (1.4) 

The author intends to show in the present paper (in the classical 
example of the strut) that also stability investigations are best handled 
in connection with the single main equation (l.2), wherein terms must be 
retained of higher. order in the deformations logically only in the ex- 
pression for the energy cf deformation. This procedure is essential 
from the practical standpoint, if the relationships are tc be investigated, 
in the poetcritical region, and desirable from the systematic standpoint, 
because it becomes clear in this manner that no additional principles are 
required. In particular, this consideration will clear up the apparent 
contradiction between equations (l.3) and (1.4). 

The calculation itself is carried out in the following manner: 

First, the expression for the energy of deformation A., is set up, then 

the differential equations for the two components of displacement are 
derived from the condition 


5 II = 0 
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and tiie question of the stability of the equilibrium position is answered 
by the restricted condition 

XT = minimum (1.5) 


Then the same problem is treated with the help of a Ritz procedure ; in 
this way the result of the stability consideration is brought out in an 
especially elegant manner.. After a thorough discussion of the usual sta- 
bility theory, it is shown in conclusion how the same considerations can 
serve for the treatment of the snap-action problem of a slightly curved 
beam. 


ENERGY OF DEFORMATION 


If the customary assumptions of the beam theory are retained — that 
for small deflections of the beam the work of stretching and work of 
bending are independent of each other and that the part of the work re- 
sulting from the shear forces is small compared with the other two parts - 
then the energy of deformation can easily be given. ■ 

As a result of the assumption of small displacements — without at first 
saying anything about the sizes of the displacements u and w (fig. l) 
relative to each other — the strains Uy. and w x (or w x 2 ) may be neglected 

in comparison with unity; that is, in a development of both quantities In 
powers of u x and w x only the lowest power need be retained. If, there- 
fore, the square of a line element of the beam centerline before deforma- 
tion is 

dx 2 

and after deformation 

dx 2 (2.1) 


(the subscripts on u and w indicate derivatives with respect to x), 
then the strain of the beam centerline (reference 1, p. 57) is 



f V 

Si 

( 

, 1 + v 

! 

+ 
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From Hooke’s law the corresponding stress Is 



and therefore the energy of stretching is 




(2.3) 


(2.4) 


The incremental strain e x due to "bending is, according to the 
assumption that normals are . preserved, ■ . 


e x ~ ~~ Zy, ' 7 XX (2.5) 

the bending energy is therefore given by 

A = | fff ( zw xx) dx dy dz - ~ ^ w xx 2dx (2- 6 ) 

o 

If eZ_= u(0) — u( 1 ) = — u(Z) Is the distance of approach of the ends, 
p = pF = - d x F the compressive force, then — (pF)(el) = y is the po- 
tential of the external forces; the total potential II (measured from 
the stress— fares state as the zero position just as the displacements 
u, w, € ! ) is therefore 


n = 


EF 


2gfV 

2 


u x + -if J a* + 


~ J ^xx 2 dx - pKZ (2.7) 


The potential per unit length - after division by EF - can be written 


Z- 


U = -S = -i 

EFZ 21 



% i S-Jta + i_ to - | E (2.8) 


1 
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From equation (2.8) end the- condition ( 1 . 5 ) 


IX * minimum (2.9) 


there is obtained all information about the behavior of the strut at and 
beyond the stability limit, _ . 


TEE DIFFERENTIAL EQUATIONS FOR. THE DISPLACEMENTS u, w 


Consider a rod the left end of which is (x. * 0 ) is supported and 
the right end gives • • (x *= l) is freely movable iri' a horizontal direc- 
tion under a centrally placed compressive force P. (See fig. I.) As 
given (that is, as the independent variable of the problem) take either 
the horizontal displacement or the compressive 

of the right-hand end stress 




According to the principle of virtual displacements, the displacements 
are to be varied under a constant load in a manner compatible with the geo- 
metrical conditions. If at a boundary point the displacement (in the pres-' 
ent problem, for example, e ) is prescribed, the point is held fixed 
during the variation, so that the work of the (unknown!) external forces 
does not remain in the calculation; if on the other hand, the force is 
given, then the (not fixed geometrically) end point is varied, and the 
work of the external load (in the present problem P 6 u) enters into the 
calculation. . . • •• 

If In equation (2.8) the displacement u is varied (that is, if the 
elements of the rod are given a virtual displacement in the axial direc- 
tion, while the boundary point or correspondingly the load is held fixed) 
there follows from the rules of the calculus of variations 


a 

8 x 


+ 


wx 


together with the boundary conditions 


u(0) = O' 
u( l ) = — e l 


u(o) « 
”x + 


0 

Wx 2 


p 

E 


x=Z 


= 0 


(3.1) 
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by varying w there is obtained (independently of whether e' or p is 
considered as given): 


8 f Wx® \ ( w x 2 \ e 

Vx §x V Ux + ” ) ~ \ Ux + ~2~ ) Wxx + 1 Wxxxx * ° 


with the boundary conditions 


(3.2) 


w(o) * w(l) =5 Wyy(o) = W„„(z) a 0 


, \ijhe exact integration of the system ( 3 . 1 ), (3.2) of nonlinear (.!) 
simultaneous .equations offers no- difficulty here. 

* _ . - ■ # 3 ' .3 

From (3.1) there follows - • . 


u x + ~~7T ~ constant = - e Q 
2 


u a u($) - e 0 x 


(3.3) 


r and with the use of (3.3)> (3.2) becomes 


v xxxx + l2 w xx = 0 


(3.»+) 


The solution of this linear equation with constant coefficients is 


w * f sin kx + g cos kx + fix + g* 


where k 1 b the positive root of the quadratic equation 


2 Go 

k 2=0 

12 


For .the determination of the six constants of integration t, g, f x , g lf 
u(o), € c there are the six boundary conditions (3.l) and (3.2): 
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8 + Si = 0 k 2 g = 0 

f sin kl + g cos k l + f±l + gx = 0 
fk 2 sin kl + ek 2 cos kl = 0 


-e 0 + - = 0 (3.5) 

E 


u(0) = 0 


and 


el = u(0) - e Q l - j dx 
o 2 


■ l 
-/ 


It is found that 


and either 


fx- - gi ~'g = u(0) = 0 


f = 0, that is, w = 0 



(3.6) 


or 


f / 0 ; v = f sin kx, k = — 

’ l 

s - ±JL_ _ e* 

° zs " 


rt f 2 jtx 

u = ~ € * x -• 8 T 8in T 


(3.7) 


There are evidently two kinds of equilihrium positions: the straight 

(f = 0) and the hent (f 4 O) . The straight position is specified uniquely 
by either' p or e, the bent by € - for the amplitude f of the deflec- 
tion can be determined uniquely from the left one of equations (3.5) 


ff 2 f 2 
■ kl s 


6 


— e 

o 


€* 


(3.8) 
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not, however, by p. From the right one of equations (3*5) there follows 
rather that for f ^ 0 a completely determined "critical" value 


p = E6* = p* (3.9) 


cannot he exceeded. Therefore p is unauited for an independent 
parameter (the situation is different in the case of the corresponding 
plate problem) (reference 4, p. 12*0. From equation (3-8) it can be seen 
that f assumes real values only for e > e*. 

In figures 2 and 3 the quantities f and p are plotted against e. 
The solid lines are for the solution (3.7), the dotted lines for equation 
(3.6). 

The result (3-9) — that the load for the buckled strut does not in- 
crease beyond p* even when the critical end shortening has been con- 
siderably exceeded - is a result of the limitation to "small" deflections. 
For the present problem this limitation is not important because here it 
was only a question of seeing that as a result of the appearance of higher 
powers of w in equation (2.8) the elastic rod can assume several equi- 
librium positions - especially that the existence of a real multiplicity 
is bound up with the exceeding of a certain. "critical" strain e » e*. 


STABLE AND UNSTABLE EQUILIBRIUM 


In the "ordinary" theory of elasticity it is necessary to consider 
only the condition 


511 a 0, that is, H = extremum 


(4.0) 


■'•Reference 5, PP. 70 ff . Also the theory of the so-called exact 
differential equation of strut buckling 

Ej/p + Pw = 0 

(reference 1, p. 280) shows that at large deflections (because of the in- 
creased demand on bending energy) there is a very small increase in load. 
From the energy standpoint, to be sure, this "exact" equation is not im- 
portant, for if w x 2 . is taken as not small compared to 1 in the bend- 
ing tern (that is, the curvature l/p is used in place of w**) it is 
necessary to proceed in a corresponding manner with the stretching term 
(equation (2 . 2) ) — unless it is assumed that an incompressible strut exists 
from the start, as is done in the theory of the Euler elastica. 
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in -the determination of equilibrium states, for the supplementary con- 
dition 8 2 IX > 0 (mechanically: the stability of the equilibrium posi- 

tion) is assured there because of the linearization (reference 1, PP* 71- 
72); in the present problem the minimal' cbndition must be set up explicitly, 
since only through 

SIX a 0, 8 2 IX > 0, ‘ that is, XI «= minimum (4,1) 


can the stable equilibrium positions be distinguished from other possi- 
ble (the uns table) positions. The concept of stability is made precise 
here by the following convention. 1 

1. An equilibrium state is called stable if for every neighboring 
state the potential energy has a larger value. 

2. An equilibrium state is called labile (unstable) if there is at 
least one neighboring state for which the potential energy is smaller. 

3. A stability limit (that is, a neutral equilibrium state) is 
spoken of when there is at least one neighboring equilibrium state the 
potential energy of which is equal to but none having potential energy 
less than that of the given equilibrium state. 

Return to equations (2.8) or (2.9): 



and perform a. variation; that is, replace u by u + 5u and w by 
w + Sw; there results, after arranging in powers of 3u, Sw and 
stopping after terms of the second order, 


1 The following definition was given, 

cideirfcally to his DVL lecture. 


in substance by E. Trefftz in— 
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,-V ’AEE = .II-(u-.+ 8u, v + 8w) r-.JJ. (h^w) « 6H + ~ 5 2 XI + . 



o 


The condition that the terms of first order shall vanish leads to equa- 
tions ( 3,1) and (3.2); the question of stability is answered by the terms 
of second order. jBy inserting for u, w the values obtained from 
SH 5= 0 there results for (3.6) (straight position) 


6 


p /s \ 

5 n j 


/ 


2? j 

Al A 




dx (^.23.) 


for (3.7) (bent position) 



j i l(- s * °° B 3 fj (®v + is ( 5 w . 


dx 


V 


l l 

2nf P nr, l Pf , 

+ — J cos — 6w x SlUj. dx + - j ( J dx (4.2 2 ) 

*■ c 7 2 o 


First the stability limit will be determined. According to the defi- 
nition given above there must be at the stability limit a state for which 
the second variation vanishes but none for which it becomes negative: 

The value zero is therefore the smallest value that 6 £ H can assume at 
the limiting point. If, therefore, 5 2 H has certain continuity prop- 
erties (the existence of which is obvious on mechanical grounds), then 
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the "characteristic" value 6 s p = Q . is at the same time an analytical 
minimum, compared with neighboring values, and the associated- ("charac- 
teristic") displacement system Su, 5v Is determined from the condition 

- e(& 2 £rj = 0 ' (4.3) 

The straight position (see equation (4.2 X )) thus reaches the sta- 
bility limit when 8u, Sw satisfy the two differential equations ' 


(5u)^ — 0 

€ 

(5w) + (5w) e 0 

xxxx ± s ' 'xx 


(4.4) 


with the end conditions 


8w(0) = 8w(Z ) * Bw^O) = 8w (l) m 5u(0) = 0 


and 


5u ( Z ) a O' 


Su X (Z) = 0 


The solution of this eigenvalue problem reads in both cases 


(4.5) 


8u =0, 5w = 5f ein 


nx 


(4.6) 


The amplitude Sf 4 0 remains undetermined and from the second of equa- 
tions (4,4) there is obtained for the critical value of 

«r n2i2 

• . crit = ~YS~ .... 

that is exactly the expression = e* by which the- branch point of 

the equilibrium was characterized; stability limit and branch point coin- 
ciae. 


The investigation proves to be somewhat more difficult 
position. The two minimal conditions read 


for the bent 


l4 
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3 nf ■ nx „ \ 

+ _ cc -• 8*j y 


= 0 


(’+. 7 ) 


■A 

Sx 


1 Sw xxx + 


6 * _ 


n s f z 2 jtx 


\ 


cos" — j 5w Y — — cos — Siu, 
I s l '' x l 1 ^ 


= 0 


for the boundary condition there is retained 


5v(0) = 5w(Z) = 8w xx (0) = 8v .^(Z) - Su(0) = 0 


xx' 


and, depending on whether or not the right-hand boundary point is pre- 
scribed or not. 


8u(z) * 0 


/ nf 
y^oux + cos 



0 


(4.7*) 


It is- recognized immediately that for f = 0, hence at the beginning of 
buckling, the bar is in neutral equilibrium, for all conditions are satis- 
fied by the solution (b.6). This result is trivial. It is not so directly 
obvious, however, that also for f ^ 0 there is a variation that makes 

8 2 H a minimum; the homogeneous system of equations with homogeneous 
boundary conditions permits of a non-vanishing solution also for f / 0. 

In fact if the variation is so performed that the second boundary con- 
dition of (h-. 7* ) is satisfied there follows from the first of equations 

(4.7) 


5u x = - ~~ COS -y 5w x (4.8) 

If this iB put in the second of equations (4.7) the latter reduces to 


i ov___ + e*5w„_ = i‘ 
xxxx xx - 


Sw, 


xxxx 


+ — Sw. 
Z* 


XX 


y 


= 0 


and this equation (together with its boundary conditions) is satisfied by 

Sw = &f sin — . For 5u there is obtained from (4.8) 
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« 2 fSf / l 2nx\ 




in which Sf again represents the (not determinable by a system of homoge- 
neous equations) arbitrary factor. 

The question of the stability 1 of the equilibrium positions below and 
above the limit can now be answered. 

T. Since the straight position (see equation (4;2 i):) is stable for 
very small s, from considerations of continuity it is so for e < e'*'.. 

2. For e > e* the straight position represents an unstable equi- 
librium state, for a variation &w, namely. 


&w * 5f sin 

can be given for which 8 2 II <0. 


3. For e > e* the bent rod is againBt the variation 


8v m 5f sin 


nx 

l 


Bu = 


3T g f 5f ( 


21 * 


V 


X + 


1 

irt 


sin 


2nx 


) 


(4.9) 


in neutral equilibrium. 1 Since the variation (4,9) (and only this) makes 
a minimum (of value zero) every other variation gives, it a positive 
value. .If in some way the special variation (4.9) is prevented, then .the 
tent equilibrium position is stable. This holds especially in the -impor- 
tant case- where not the force but the displacement of the end point is 
prescribed; for the displacement system (4.9) is in fact excluded by the 
boundary, condition 6u(Z) = 0. 

At this point a result' discussed later in section Connection between 
the Ordinary Investigation of Stability and the Procedure Presented Here, 
should, be emphasized. The behavior of the rod beyond the stability limit 
is Different depending upon whether the load' or the displacement is 

;• .- This result, is naturally connected with the assumption that the 
strut adheres strictly, -to .the. law of deformation established by’ the ex- 
pression (2.7). ” 
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considered as the prescribed quantity. More note worthily, however, the 
"behavior at the stability limit is not affected by this difference. For, 
although there are the two different boundary conditions (^.5) 


8u (l) = 0 

.and 

5uj (l) c 0 

they both (together with the differential equation &U y r = 0) lead to 
the same result 

&u = 0 


That is, it makes no difference whether a motion of the end points in the x— 
direction is "permitted" or not: during the buckling they do not move. 

Thus the result is arrived at that the two mechanically entirely different 
problems: buckling under constant load and buckling under constant end 

shortening, lead to exactly the same critical state Su, Sw, e cr it^ Pcrit.* 


INTERPRETATION CF THE RESULTS WITH AID OF THE RITZ METHOD 


The results of sections The Differential Equations for the Displace- 
ments u, w and Stable and Unstable Equilibrium can be illustrated very 
elegantly if the variation problem is turned into an ordinary minimum 
problem by the use of the Ritz method. In the case of the Ritz method 
to be sure nothing certain can be said about the question of stability, 
since from the start only' quite definite displacements are considered and 
therefore nothing general can be concluded about the sign of the second varia- 
tion; nevertheless, with a Judiciously chosen deflection system the 
question can be answered with great probability or the answer made very 
plausible. In the present, especially simple case -the earlier results 
are found again exactly. ' • 

As a set of displacements - satisfying all boundary conditions are 
chosen, the solutions of the differential equations (3.2) and (3.l) 


x This is a peculiarity of the problem. In general, the critical load 
depends upon the. boundary conditions — whether are prescribed forces or 
displacements, the system is supported, or guided, or built in, etc.; for 
the "minimum" - variation, from which the critical load follows, differs 
according to the boundary conditions prescribed by the data of the problem. 
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. wy 2 _ . « 2 f 2 

v mf 8in T* "x + = ~ e + Tzi 


(5.1) 


Then the minimal condition 

l 


TL*±f U + dx + g fv^dx-f 

21 J \ 2 / 2 1 J x* E 


«V « e f a 


_ , (e ~ e*) + — . — — e = minimum 

321 If] 2 E 


(5.2) 


furnishes an equation for the "free value" f as a function of € 


is a t£LL* _ s + ££. ) „ 0 

Sf 2Z 1 2 v hi 3 


(5-3) 


The relationship between load and end displacement is obtained by means 
of the etre s e-strain equation (2.3) from the second of equations (5.1); 
it becomes 



n 2 f a 




(5.4) 


There are again the two possible solutions 


f « 0, p = Ee (5.5) 

and 

' jr 2 f 2 ' 

f f 0, hence — 3 — - e — e*, p = Ee* (5.6) 

hi 

The question of stability is answered (with the above-mentioned 
limitation) by the second variation. As hitherto, two cases are distin- 
guished: 

1. If the force is prescribed, € therefore left open, there are 

two displacements to be varied, and the sign of the expression 
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p A 

5 II = 


5 g n 

df 2 


(5f) : 


'.+ 2 


5 2 


IT 


<5f3e 


SifSe + 


5U 

de 2 


(Se) : 


(5.7) 


must be investigated. This expression is a quadratic form in Sf , 8e 
with the coefficients 


a g n 

df 2 



+ 3 


2 2 
W f 

4l 2 


) 


8 g n _ fir a 

Bfde ~ 2l 2 


Since 


n 2 ~ 

— -S. > 0 
& 2 u ' 


P-B. = 1 

de 2 


thiB expression is positive definite (that is. 


never negative) as long as the discriminant 


A = 


00 _ o 

d II d^II 
Sf 2 de 2 


a g nV 

dfdey 


2Z C 


* _ 


€ + 


2 2 
Jl A 


4l £ 


(5-8) 


is greater than zero. This is certainly the case below the critical point 
(c < e*); here the system is therefore stable. On the other hand, A< 0 
for f = 0 and e > e*; that is, the straight position is unstable above 
the critical point. Finally, the bent position is neutral, because by 
(5.6) the discriminant vanishes for f ^ 0, This result agrees with that 
of the previous section and is an expression of the fact that a buckled 
etrut in the elastic range can be bent arbitrarily farther without in- 
creasing the load, 

2. If the end displacement e is prescribed, only the quantity f 
need be varied and it is found that: ' 


8f 2 2l a \. 4z 2 / 


(Sf ) 2 . 


(5.9) 
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From this relation it follows immediately that 

for e < e* > 0 

- for e > e* and, . f ^ 0 S 2 ^ > 0 . .. 
for e > e* and f '• 0 5 2 H < 0 * , 


that is, the straight position is stable for e < e*, unstable for 
c > e*, the bent position, as soon as it is mechanically possible (hence 
for 6 > e *), always stable. Figure k shows the energy relationships in 
this second case (prescribed displacement e = as* of the right-hand end 
of the rod). Plotted as 


or 


tk. . i e 4 

e* s 2 



- 1 ) + 


a 


2 


as a function of f(^or | - &-J with e for a « as a parameter. 

(The potential of the external forces is not included, because it is not 
affected by the minimal condition with respect to f.) It is seen that 
the straight position (f =0) is an equilibrium position under all 
circumstances, for all curves start with a horizontal tangent. • For 
a < 1 associated with f — ■ 0 is a minimum, for a > 1 a maximum;- the 
curves for a > 1 have further to the right also a minimum, whereby the 
bent position f ± 0 is characterized as a (stable) equilibrium position. 
This figure shows especially well the "type" change of the curves in the 
transition from the sub-critical to the supercritical region: the co in— 
ciding of maximum and minimum for a =* 1 . It is also clear here that, 
although at the moment of transition the displacements are small , the 
behavior of the body at the stability limit is nevertheless determined 
by the "possibility" of greater deflections, expressed mathematically, 
by- the existence of the f— terms of higher order in the expression for 

V 
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C DIRECTION BETWEEN THE CSEBUfflY INVESTIGATION OF STABILITY 
AND THE PBOCEDOHE PRESENTED HERE 

The ordinary stability' theory Is limited to an investigation of the 
critical point. It was seen that the critical point is characterized by 
two energy conditions. The condition 


• : 5H = 0 

for any variation Su, &v, Sw . (6.1) 

characterizes it as an equilibrium position in general; the condition 

sfn = 0 

for a characteristic variation Su, Sv, 8w (6.2) 

as the critical one. Or, the critical point is distinguished by the fact 
that there a variation of the state of deformation can be made for which 
the potential H remains unchanged to terms of the second order. Now 
in practical buckling problems it is usually a question of the transition 
from a very simple (often Independent of the coordinates) initial state 
of stress to a comparatively very complicated one. It Is therefore custom- 
ary to specify the initial state 'of stress directly without recourse to the 
definition in terms of energy (6.1) and to proceed with the variation im- 
mediately in regard to the determination of the- second state. There then 
remains as the single important condition the statement (6,1), which can 
be expressed in the form of a method-of-procedure as follows, for example i 
Consider a. system of infinitesimal distortions superimposed upon the 
critical state of deformation, collect the parts of the potential . energy 
II quadratic in the added displacements and set, the sum equal to zero. 

That .such a procedure is at all possible rests Upon the fact that as a 
result, of the "large" initial stresses two types of quadratic term arise: 
an (always positive) part that represents the work done by the stresses 
caused by the added displacements, and a second part that comes from the 
work done by the stresses already present upon the quadratic part of the 
added displacements. (See reference 2.) 1 

The fact that this statement--of -procedure concerning the vanishing 
of the quadratic members is nothing more than the extended principle of 
virtual displacements ("extended" in the sense of the statement about 
8 2 U) does not come out clearly in the applications mostly for three 
reasons . 

equation ( 2j. ) , for example, the last two terms represent the 
first type and the first term, tbo second type. 
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1. Since a confusion with the very simple initial state is in general 
not to he feared, it is possible' to dispense with the designation Su, 

8v, 5w and write more briefly u, v, w for the added displacements. 

This manner of writing does. not express the fact that the added displace- 
ments are to be not only small in the sense of the general hypotheses of 
the theory of elasticity but also infinitesimal in the sense of the calcu- 
lus of variations. - ... 

2. In close connection with the above, in considering the energy it 
is customary to start not with the total potential II but directly wiih 
the energy changes (appearing as the, re suit of u, v, w) and to designate 
these changes by 1 A, _V instead of by 5A, 8V; the (extended) principle 
,of virtual displacements becomes in this manner of writing 


Aj + V = 0, or even A^ *= A a (6.3) 


since the potential difference —X also represents the work of the 
external forces on the infinitesimal 2 displacements u, v, w. Equations 
(6.3) can be put into words as follows: For the virtual displacement u, 

v, w, through which the original equilibrium configuration goes over 
into the neighboring ("buckled") configuration at the critical point the 
internal energy A^ taken up by the system is equal to the work done by 

the external forces A a taking into account the terms linear and quadratic 
in u, v, w. By this formulation the two conditions (6.1) and (6.2) are 
combined into one; a procedure in which there is the danger of losing 
sight of the difference between the (holding for any equilibrium position) 
principle (6.1). and the (characterizing the critical position) extension 
( 6 . 2 ). 

3. As the, proper equation for the determination of the critical 
system of virtual displacements u, v, w there follows (see sec. Stable 
and Unstable Equilibrium) from (6.2) and the added requirement 

S 2 H > 0' for all other- 8u, ■ Sv, 5w 

the condition 

. 8(5 2 H) = 0 (6.4) 


To distinguish them from those used earlier, the quantities usually 
written A ± , A & , V are designated by Af, A a , V. 

2 The second form of the law (6.3) therefore does not represent the 
special energy law A i - V by which is expressed the fact that for 

conservative systems the external work introduced hy the transition from 
the initial to the (not neighboring!) final state is stored up as elastic 
energy in the body. 
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If it is agreed to consider in A & , V_ only, the (alone essential for 

the critical behavior) quadratic terms > the condition (6.4) is written 
in the form 


5(A i '+ y) =0 
“ A^) - 0 


(6.4*) 


.This form, which is only a natural consequence- of. the original agree- 
ment to mite u, v, w instead of. &u, bv, Sw, makes it quite clear to 
what extent the simplified manner of writing can: lead .to conceptual errors. 
For the statement (6.4*)/ aside from the deceptive formal agreement, has 
nothing to do with the principle of virtual displacements (l.2) or (1.3): 

The principle (.1.2), in content the same as the energy law (see the 
Introduction), answers the question of the equilibrium positions under 
given loads (or edge displacements), and equation (l. 3) is a special form 
of the same principle possible only in the realm of linearized elasticity 
theory besides being very inexpedient 1 ; equation (6.4'J, on the .other hand, 
in content the same as the. minimal condition (6.4) concerning the behavior 
of the quadratic terms at- the stability limit, gives the second equilibrium 
position possible at the -branch point and the sought— for value of the load 
at which the equilibrium begins to be many— valued. 

The difficulties so far discussed- were difficulties in interpretation 
arising from the symbolism of writing. There is another, more factual 
circumstance that makes the question compiex, especially difficult to see 
through. It was seen in section, Stable and Unstable Equilibrium that far 
the rod there were two independent equations (4.4), with the likewise 
independent boundary conditions (4.5), for the two added displacements 
5u, Sw (which here would have been written u, w). From them it was 
concluded that u vanished identically. This result — and correspondirgly 
u = 0, v = 0 in the case of plates — makes possible, when (as is tacitly 
done, in the stability theory of a bar) J.t is. presupposed as known, a treat- 
ment of the problems of bar and plate stability deviating from the general 
methods of stability theory depicted above. Since, however, bars and plates 
are the most well-known problems, being analytically the most tractable, 
frequently ideas that were developed there are erroneously brought into 


1 So, for instance, for the compressed strut below the critical point 
twice the external work can be. written in three forms: £s 2 , pe, p 2 /fe — 

which is to be varied (with respect to . e!)? The second form is meant, 
but as a result of writing 2A a in the place of —V that is no longer 

uniquely discernible. 
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more general stability problems. It is therefore necessary to examine 
more thoroughly the various — special — interpretations that can be given 
to the occurrences at the stability limit in the case of bars and plates. 

First of all outline the method by which it is necessary to proceed 
according to the directions formulated at the beginning of this section. 
If it is assumed a virtual, displacement u, v at the critical point, 
then, as can easily be seen 1 the strain of a fiber to terms of the second 
order is given by 

2 

w x : ' . 

€ x c ”x + — “ zw xx = e x - zw xx 


Therefore the terms of second order are: in the work done by the added 
stress 




dx dy dz 


in the work done by the already present (critical) compressive stress p 

2 


Wv 


. , P dx dy dz 

JJ 2 

(The external force does work - Pu(l); this makes no contribution of 
the second order.) After integration over y and z there results 


EF 

2 


l l 

^ y\. 2 . ,2 r_ 2 


Ai = -7T l / u** dx + l 2 J w ^ 2 dx - | J w 2 


000 

and the conditions (6,3) and (6.4 f ) become 


dx j A a = . 0 


(6.5) 


r to ™. 2 dx + i 2 r v. 

■£ -L - 


2 5 

'XX E J ’x 

'O .0 


is 


dx = minimum = 0 


(6.5*) 


The term u x /2 goes put in the expansion of the radical (see . 


equation (2.2)), and the expansion of the curvature — 

P 

would give terms of the third order. 


w xx 




y 
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The expression coincides perfectly with the earlier expression (b.2j.); 
therefore the same differential equations and boundary conditions, and es- 
pecially the result U 5’ 0 are obtained entirely independently of whether 
.' a motion in the x-^lirection of the right-hand end point during the buckling 
is permitted or prevented. This double result (that u * 0, and that the 
. buckling is independent of the condition u(l) = 0 or u^(l) = 0) makes 

possible the two following "customary" interpretations of the buckling 
process-. 

The first procedure consists In considering Instead of the "natural" 
problem, buckling under fixed load, the problem of buckling under fixed 
end point and at the same time (what seems almost a natural consequence 
of this stipulation) assuming from the start the vanishing of u also in 
the interior. Hence there is supsr imposed upon the straight position 
w = 0 a purely transverse displacement as a variation, keeping in mind 
the presence of the still unknown longitudinal compressive stress ~Ee 0 . 

According to equation (2.2), as a consequence of the change in length 
connected with the transverse displacement, the following stretching 
energy is released 


l l 

P , _ EFe 0 r g 

Ai = / (~EFe 0 ) e x dx = - -5— / w x dx 

v J o o ■ 

at the same time a bending energy 

l 

EJ P g 

• . , ^-i ~ 2 ■ J w xx ^ 

must be added. This interplay between the two types of energy X and hence 
the two equilibrium. positions) takes place when the values of Xj_ and 

are numerically exactly equal; that is, when 

l ' l : ■' 

Ai (i 2 j ^ w/ dx J » 0 (6.6) 


From the additional condition SAi =0 there follows as above the sine 
equation for ,w. This procedure thus leads to the correct end result 
. without, however^- permitting, a" guarantee.; of really having 1 found the mini- 
mal buckling load. For the "restraint" assumption u s 0 limits the 
number of possible variations, and that it leads to the correct buckling 
load for the rod (and plate) requires at- least a supplementary verifi- 
cation. 
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More Important — 'because in a still more special way a peculiarity 
of the rod. and plate — is another manner of thinking, which is almost 
universally made the basis of derivation of the buckling equations. With 
reference to the natural huckling process, the "boundaries are considered as 
movable; however — and this Is the characteristic mark of this method — 
they cannot he allowed virtual displacement u (or , u, v — which, as haB 
been observed, would subsequently become zero) but are given a displace- 
ment that is of a higher order of smallness (compared with w). 

In the case of the beam it is customary to start this procedure with 
the assumption that no additional stretching energy is taken up during 
bending; that is, that the bent beam has the same length as the straight 
one; it follows therefrom that as a result of the bending the ends must 
approach each other by an amount 

l 

| yW 2 (6.7) 


(which in fact is of the second order in wl), so that the external forces 

l 

p F P 

do the work - p Q Fu j / w x s dx. Now by formulating the equality 

2 J o 

of inner and outer work (wherein by inner work is to be understood only 
the bending energy) 


= A & or 


EJ 


w xx dx ~ 


EFe, 


w x 2 dx = 0 


and assuming as above the minimal property of this expression, this pro- 
cedure leads to equation (6.5*) - naturally likewise without the u-term. 1 


In the assumption-. (6, 7) there is an inconsistency: It cannot be 

assumed a priori that a strut that changes its length elastically below 
the buckling limit suddenly ceases to do so beyond it. (in reality it 
changes its length by quantities of higher order . ) It is more logical 
to consider a perfectly incompressible - rigid against extension but 
elastic in bending - strut, for which the two hitherto independent dis- 
placements u and w are related from the beginning by the (geometrical ) 
assumption 

W-jr 2 

+ ■— » 0 ( 6 . 7 *) 

(Continued on p. 26) 
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Since the procedure of equating the stretching energy to the external 
•work cannot he used in the case of the plate, a special auxiliary idea 
has been used there in order to preserve the conceptually so similar idea, 
that the boundaries are to move. (See reference 6.) 

Without connecting the displacements u, v, w with each other nu- 
merically proceed, in this method, from the assumption that u, v are 
of a higher order of smallness than w; that is, consider u, v not 
as really .independent virtual displacements but as connected with the 
transverse displacement w by the order- of- magnitude condition 

u x * w x 2 (9*0.) 

.According to equation (2.2) there is obtained for the work done by the 
critical stresses cy, ay, r on the displacements u, v, w to terms of 

the "second order" . . 

rr/ v x 2 \ rr / v x a \ 

Ax = ^X s j I yu x + — ) dx dy + o yS j j (^v y + — j dx dy 


( 6 . 8 X ) 


(6.82) 

(Continued from p. 25) 

Such a strut permits no deformation at all below the critical, load; above 
it takes on only bending energy, which is furnished by the exter nal work 

1 ■ r l 

2 ScF J 0 v x 2 a* 

Since up to- the critical load no elastic deformation at all has taken place, 
the two laws 

= 0 

are here in content completely identical. 


Ap + I = E'CAj + y) 


2. 

5 (A- 


S (A, + V) 


rr/ \ 

+ TB 11 V + v x + w x w y J ^ 

the bending energy is, as always, given by 

4i = lo(l ~~)ff (AW)S “ 2(l_v) ( w xx w yy ~ v xy S )_ dx 
(s = thickness) 
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•Both parts together must he equal to the work A a ..of the external force 

in the sense of equation (6.3).' The external woirk can' now (and this is 
the essence of Rdissner 1 s idea) he ©xpi?ee8ed generally -in a wry simple 
manner, if it is femeribesed that the straight position is an equilibrium 
position and that therefore in every* virtual' displacement A a - A^. On 

taking the special displacement u* = u, v* ■ v, v* = 0 (with the bound- 
ary displacements, u^* * .u^,. v R * = v^), A & = A^*; therefore 


Ct X a X + 0 y v y + T (^ u y + v xJ j dx dy 05. 83 ) 

and now on collecting terms in (6.8), the u- and v-terms cancel out; 
there results the well— known Bryan plate equation (reference 1, p. 293), 
exactly in the form obtained also under the assumption of purely trans- 
verse displacements and immovable boundaries. 

The advantage of this method is that it offers the possibility of 
. formulating- exactly, the related presentation of a solution of the buckling 
process by a. boundary displacement. Its disadvantage is a double one: 

The emphasizing of the boundary, displacements gives the Impression that 
the participation of the external work is. universally important in a 
buckling process, which, as has been seen, is not so. But besides this it 
is important' for the entire consideration. Just as for that of Bryan, that 
u, v are of the second order with respect to w, which must be known 
somehow beforehand; 1 therefore the interpretation of the external work 
a x u x • * "* 8,13 a contribution of the second order, is not transferable to 

more general buckling problems. (See reference 6.) 

To summarize briefly the result of this section: In considering 

the critical point it is customary to dispense with the correct method 
of writing the virtual displacements 6u, 5v, 5v in favor of the more 
convenient u, v, w; thereby the connection between the customary 
stability criterion and the principle of virtual displacements is con- 
cealed. To be added is that the stability problem of the rod and of the 
plate permite a special -treatment which rests' upon the fact that at the 
critical point the tangential displacements u, v and the normal 


1 The very obvious conclusion, that can Just be seen from the form 
(2.2) of the strain that u x end w x 2 must be of the same order, is 

not tenable; for an equation of the type (2.2) holds, for instance, also 
•for the longitudinal fibers of a cylinder, and yet here and w x 

can become ' comparable because the tangential displacement v, which is 
of the same order as w, is linearly coupled with u through the shear 
and the transverse contraction. 
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displacement w are of different orders of magnitude. Since, however, 
the rod or plate problem, as the' analytically simplest, is at ’the same 
time the best known, the need easily arises of transferring methods Of 
thinking successful in these problems to more complicated problems, which, 
as was to be shown, is not possible. 

THE DUBCESCHLAG - FBC8LEM OF THE SLIGHTLY CURVED BEAM 


In section Stable and Unstable Equilibrium, it was shown that for 
the Euler strut the instability point (defined by 5(S 2 H) = 0) coincided 
with the branch point of the equilibrium. Branching problems are, however 
not the only kind of stability problem; a second class, which is Just as 
suited to the energy definition, of stability as are the’, branching 
problems, comprises the so-called Durchschlag problems. 

In the Durchschlag . problem the critical load is designated as that 
load under which an (infinitesimal) displacement of the point of applica- 
tion of the load is possible without an increase in the load, for which — 
as in the branching problem - there are therefore two (infinitesimally 
close) equilibrium positions. Above the critical point an increasing 
displacement is in general accompanied, by a decreasing load — the state 
is unstable, the system "snaps" into 1 or falls into a stable configuration 
Prerequisite for such a phenomenon is a nonlinear relationship between 
force and displacement even in the stable region. 

The simplest Durchschlag problem is that of a slightly curved 
beam under a transverse load. (See reference 7, ) 

If the ends of the initially curved beam are prevented from dis- 
placing (fig.' 5), then connected with the deflection caused by the 
transverse force Q is' a shortening of the axis of the arc, as the 
result of which a horizontal force H is made to act. Because the 
effect of this (very large) congress ive force upon the equilibrium of 
forces in an element of the beam cannot be neglected, there arise 
phenomena related 1 to the buckling process in the Euler strut, insta- 
bility phenomena. 

Without carrying out all the details of the calculation (presented 
completely elsewhere, reference 8) the principal method of solution for 
this stability problem will be briefly outlined. 


1 Xn a manner, similar to that- in which a strut compressed beyond 

the Euler- limit at the least disturbance snaps or falls into the bent . 
position. 


* 
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By a process that follows very closely that -carried out in section 
Energy of Deformation, is obtained, with the notation of figure 5 (w 
now taken positive downward), for the potential energy 


n = 


EF 


/ [<S 1 


„ V . .a 2 


+ + 1 w - 


XX 


dz - Qf ; (7.1) 


From 


■ • : . r a 

6IX 2 EF jf + ~~ - W^w Q 6u x dx 


+ EF 


f _(% + . — W x w x)( w x - W x) 6w x + l2w xx 8v xx 


dx - QBf 


there are obtained the two equilibrium conditions expressed in terms of 
the displacements: 


S / v x 2 

+ T -V* 


or. Integrated once. 


( 


"x- 


w x 


W x w x 


= constant = — h 


(7.2) 


and, with the use of equation (7.2) 


1 W XXXX + ^ W XX “ 


(7.3) 


also the boundary conditions . 

u(0) = u(i) = w(0) * w(Z) = Vxx (0) = Vxx (l) = 0 

w xxx( J /2 + 0) - w ^(2/2 - 0) * Q/EJ 


(7.^) 
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which together with the continuity requirements . for 


u, u x , v, v x , w^ 


at< x = 1 / 2 give the 12 conditions that are necessary for the evaluation 
of the 2x6 constants of integration in the 2 regions x ^ l/ 2. 

The physical significance of the constant h can "be recognized di- 
rectly from (7.2): On the left-hand side is the stretching of the middle 

line of the arch; therefore, to a factor l/EF h is equal to the hori- 
zontal force H, and (7-2) expresses the. equilibrium condition that H 
does not vary with x. 

Just as in the case of the Euler strut the constant of integration 
h = H/EF of the first equation enters into the second as the coefficient 
of the unknown w; that is, the system (7.2) and (7-3) is nonlinear, 

(See equation ( 3 .^) or ( 3 . 2)).. Nevertheless, Just as before the exact 
solution can be given in this simple case without difficulty in terms of 
the at first unknown 


l 2 h H 


a = 


2.2 

n i 


= — , - H* « 

H" 


EJit 2 


JtX 


for example, for W = f sin — r 


a 


w = f, 


O 2 
a 


JtX 

sin — + 


0 , 1 ' 


airx 

sin — 5 — 

1 curx 


l 2n 2 a 2 EJ cos ~ 


(x < 1/2) 


(7.5) 


and from (7.2) by another integration taking into account the boundary 
conditions u(0) = u(2 ) = 0 there is obtained subsequently a transcen- 
dental equation for the dependence of a upon Q and f 0 of the form: 



;(7.6) 
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The da termination of the critical load can he carried out in two 
basically different ways. 


The first, method (see C . B , Biezeno, reference 7, pp. 21 ff) proceeds 
from the condition bQ/bf *< 0, herein the relation between Q and f 
is established by equations (7.6) and (7.5) for x « 1/2 


w(z/2) * f - f c 




2jr 2 a®EJ 



(7. ‘5*) 


A second method proceeds by way of the energy criterion (^.3) , 
For & 2 XI is obtained from equation (7.1) after writing for brevity 3- 


- »X^)* (»x - -x) 8 •] • - [ ~ - ( »x : ~ "x) 


= — w (x,a) (7.7) 


the expression 
l 


8 2 H «EF. J' . jjsuj.) + 2 (v x ~ W x ^6u x 6w :jc -u(x,a)(^6v x ) + ix 

- J • (7.3) 


The particular displacement system 6u, Sw by which 5 2 n is Just made 
equal to- aer© is obtained from 


6(6 2 H) * 0 


that is, from the two homogeneous differential equations 


. ’SUtt- + — 

** dx 


SWt, 


and 


~ <ic ( Wx “ w ^ 5u > 


(j T x w x^) 

i^wxxx + w(a,x)5w x J * 


dx 


• (7.9) 


0 


■ Therein v is at first according to (7.5) a function of the two 
parameters a 'and Q. Q for example, is considered as eliminated with 
•the help of equations (7.6) and (7.5*). 
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with the homogeneous boundary conditions . . 


5u(0) * 5u(l) = 5w(0) = 5v(i) * Sw^Co) » Swyy ( 1 ) = 0 

* • . ' < 

The desired critical a-value is the lowest eigenvalue of the equations 

(7.9). 


APPROXIMATE DETERMINATION OF THE SNAP LOAD 

Because of the great- mathematical difficulties that equations (7.9) 
-present, the second method outlined is not suitable for an exact treat- 
ment of the problem but is well suited - and therefore that procedure 
will be considered here - to an approximate treatment by the method of 
Eitz or Galerkin. 

This procedure can be started at either of two points: either, 
make a Eitz approximation for 8w in (7.9), determine the correspond- 
ing Su from the first of equations (7.9), and following Galerkin from 
the condition 8(5 II) * 0 obtain a (transcendental) equation for the 
determination of a; or - very much more simply, if : also necessarily 
with a - corresponding loss in accuracy - introduce at the start 
a Eitz approximation for w itself in place of (7.5) into the expression 
(7.1). 


It is well to use the second method but only indicate (reference 8) 
the course «f the calculation. If again 


W * 


sin 


jtx 

T 


is chosen and as a Eitz expression 


w = f x sin — + f 2 sin — (8.1) 

then all boundary conditions are fulfilled, except for the one dis- 
continuity-requirement '(7.4), the violation of which is. however, un- 
important. Further, by satisfying exactly equation (7.2) (obtained by 
variation with respect to u) and calculating the horizontal force 
H. from (7.6), the integral in (7.1) can be evaluated and II is 
.obtained as a function- of the amplitudes f Q , f x , f 2 or the dimension- 
less parameters " 
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£ fx 

- “S A x = — , A e =» — (l * radius of gyration) 

--■ - - i i i 


in the simple form 

~ 42 4 n s , s i / V 2 9\ 2 

^ = it 4 i 4 EF2 “ ^ + -^Ag + ~ ^ ~ **7; 2Ag j} - qA 0 Ax (8.2) 


where 


_ 1_ 4l 3 Q_ _ 42 Q2 2 _ 42 _Q_ 

V it 4 ! 3 E3P n e fp EJjt s ” jr 2 f 0 H* 


( 8 . 2 *) 


The equilibrium equations read 


dll 


= 2A a + (A 0 — A x ) 


(* 0*1 


A, 2 


- 2A 2 ) - qA Q = 0 


?)TT /■ A, 2 

~ m 3 2A a - 4A 2 ^A 0 *x - - 2A/^= 0 


(8.3) 


They are in the two unknowns A x and A 2 and of the third degree; 
nevertheless a complete discussion is possible without numerical calcu- 
lation, because the second equation may be written as a product 


^=(*0^1 “ ~7T' ~ 2 ^a 2 - 8y « 0 
Therefore the cases can be distinguished 

*e m 0 

2A X + (A q - Ax) (* 0 Ax - J = qA c 
and A 2 5/ 0, that is, . 

2*2 = A 0 *x — . A x /2 -r 8 


(8.4) 


8A C — 6Ax = qA 0 


(8.5) 
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the discussion of which no longer requires .any labor-/' The first system 
provides a symmetrical deformation,- the second a superposition of a sym- 
metrical and an antisymmetrical deformation. The corresponding horizontal 
compressive forces become 


Hj = H* 


^ 0^1 — ^r g / 2 
2 


(8.4*) 


H- 


II 


A 0 Ax “ *i 2 /2 - 24 2 2 
= H* : = 4E* 


(8.5*) 


The critical load q^.^ is found from the condition 


2 * 

S IT = 


-.2 ~ 

a n 


(6\ f + 2 


rn 


(S^x &4 2 ) + 


a 2 n 


( 5 ^ 2 ) = 0 (8,6) 


As long as 


V 9 A 

aAx 2 


and the discriminant 


a 2 n a 2 fr _ / a 2 n f 

8Ax 2 SA 2 2 \ aAx8A 2 / 


are greater than zero, 5 2 n as a positive-definite quadratic form in 
SAx and &A 2 cannot be zero for any combination of these two variables. 
Vanishing of the discriminant characterizes the pair of values A x , *s 
for which there is exactly one combination 6A X , SA 2 for which 5 2 II 
becomes zero but none for which it is less than zero. The condition 


^ 2 N p a 

an an 


aAx 2 aA 2 2 


a s n 

3Ax3Xs 


2.. 


= 0 


(8.7) 


gives, therefore, the stability limit and together with the two equations 
(8.3) determines the three unknowns A x , and q^-it* case 

(8.7) reads 
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+ (* 0 - *i) 2 - X x - ~-j 

+ 2^2^ J — l6^2 2 (\> ’ ~ ) 2 ** 0 


which condition becomes 


with 




8 - (\A - )] [ 2 + (*© - * x f - ( \ Q \ - 


= 0 


and with 


V 


2* 2 2 . Vl “ -§- - 8 


Wt Vi - V/2 - 8 ) = 0 


There are’ therefore two sets of values for V V and q crit 

■Ks = 0, = 7v 0 ~VV - V Icrlt = 2 + ^ ^^) 3/2 (8.8) 


*2 ® °» \ = *0 - \/^o 2 “ 16, = 2 - -A 0 2 - 16 


(8.9) 


hy which a critical state of the elastic system is characterized. 

The physical significance of the relations ( 8 . 8 ), ( 8 . 9 ) and espe- 
cially the -(in this case unstable) behavior above the critical load will 
not be pursued in detail (see reference 8 ). (See fig. 6 .) There the 
load Q is plotted against the deflection f x of the point of application. 
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with the initial amplitude f 0 as a parameter. For f^/i « - 2 Q 

increases monotonically with fj.; instability is not possible.,. In the 
region 

2 < \ ^ /22 ( 8 . 8 ») 


under the critical load given by (8.8) there enters definite instability, 
increasing deflection f* without increase of load. Accordingly the 
beam snaps - under constant load - until it finds a stable configura- 
tion 1 at Since > )\ 0 the beam is now convex downward; it also 

can be seen clearly that the system now must be stable with respect to 
an increase in load; a further deflection results in a longitudinal 
pull. (See equation (8.4*).) 

For 

* 0 W22 N (8.9*) 

the critical load is given by ( 8 . 9 ). Before the external load can 
assume the value (8.8), the longitudinal compression according to (8,5*) 
reaches the value 4 h*, that is, the second Euler load, under which the 
strut assumes the S-shape configuration 4 0* It snaps again into 
a stable position 7\xS this time, however, passing through an un~ 
symmetrical deformation. At the critical load there appears a branch- 
ing of the elastic equilibrium; figure 6 shows the. two branches of the 
Q - fx curve, both of which however - end this is the noteworthy differ- 
ence from the Euler problem — are unstable. 

For further details see' the publication referred to. Here it was 
just a question of presenting the chain of ideas that led to the deter- 
mination of the critical loads (8.8) and (8.9), in order to show the 
application of the general stability criterion (k,3) to a stability 
problem of an- entirely different kind. 


1 That is, at first vibrates about Ax’ as a stable equilibrium 
position. 
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TRANSLATOR'S NOTES 

Equation 3*7, last part 

Trane, notei It appears that this equation should he 

* f2 , 2nx ^ 2 ' 

u “ - 6 x ~ C — 8ln — ■ — x * 
■ 8 i i 4: 

Equations (7*5) and (7.5*) 

. • Ql 3 

Trans, note: It aopears that the term — r^r~ ~ 

• 2it 2 a S! EJ 


should he 


Q* 


2it 3 a 3 EJ 


Equation (8.8) 

Trane, note: It appears that 


2 + - — i — — ) ' should he 

V 3 / 



Page 24 
Trans, note: 


A, 


[«. 

j 

o 

f v ™ 


E J / 2 

It appears that Aj_ = -g- / dx should he 
, o 

7, 

EJ 


Page 36 

Trans, note: It appears that (8.4*) should he (8.5*) 
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Figure 2,~ Amplitude f against e. 



Figure 3.- Load P against e. 
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Figure 4.- Variation of energy of deformation 
v/ith the amplitude f = 2i£. Edge 
compression e = ae* as a parameter. 


Moment of inertia J , Section F 

Radius of gyration i, Elasticity modulus E 



Figure 5.- Slightly curved team under 
transverse load Q. 











